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Abstract. We present results from three-dimensional ideal magnetohydrodynamic 
simulations of low (3 compact toroid (CT) injection into a hot strongly magnetized 
plasma, with the aim of providing insight into CT fueling of a tokamak with parameters 
relevant for ITER (International Thermonuclear Experimental Reactor). A regime is 
identified in terms of CT injection speed and CT-to-background magnetic field ratio 
that appears promising for precise core fueling. Shock-dominated regimes, which are 
probably unfavorable for tokamak fueling, are also identified. The CT penetration 
depth is proportional to the CT injection speed and density. The entire CT evolution 
can be divided into three stages: (1) initial penetration, (2) compression in the direction 
of propagation, and reconnection with the background magnetic field, and (3) coming 
to rest and spreading in the direction perpendicular to injection. Tilting of the CT is 
not observed due to the fast transit time of the CT across the background plasma. 
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1. Introduction 

It is important to deliver fuel into the core of a tokamak fusion plasma to maintain 
steady-state operation, achieve more efficient utilization of deuterium-tritium fuel, 
and optimize the energy confinement time [l]. Several fueling schemes have been 
proposed, such as edge gas puffing, pellet injection p|, and compact toroid (CT) 
injection [3]. Among them, CT fueling is considered to be the most promising method 
for core fueling because the injection speed via this method is far higher than those 
of the other methods. Although extensive worldwide efforts have been devoted to 
study CT fueling theoretically [U, El H], numerically [5], El [TJ E], and experimentally 
[9l[IOl[IIl[l2l[T3l[Tl[T5l[l6l[I71[T8l[l9lE^ dynamics of 

core CT fueling of large devices like ITER (International Thermonuclear Experimental 
Reactor) [2H] is not well understood. CT injection has the potential to deposit fuel in 
a controlled manner at any point in the machine, from the edge to the core. Tangential 
(toroidal) injection can impart momentum for improving plasma (3 and stability p9]. In 
a burning plasma device with only radio-frequency (rf) for auxiliary current drive, a CT 
injection system may be the only internal profile control tool for optimizing bootstrap 
current and maintaining optimized fusion burn conditions. CT fueling also provides a 
good chance to study core transport in present machines, helium ash removal, and Edge 
Localized Mode (ELM) [30] control. 

In this work, we employ a simple idealized model of a low (3 CT propagating into 
a uniform slab plasma with a uniform magnetic field perpendicular to the CT injection 
direction, mimicking CT fueling into a tokamak with infinite aspect ratio. This model 
helps us identify different regimes of operation in terms of CT injection speed, density, 
and magnetic field strength, as well as understand the essential physics occurring during 
CT injection. More realistic scenarios, including the use of realistic tokamak profiles 
and geometry in the background plasma, as well as high /3 CT's and dense plasma jets, 
are planned for follow-on research. Compared to past work on CT injection simulations, 
we have investigated new regimes especially in terms of higher injection velocity and 
a more ITER-relevant ratio (at least for low f3 CT's such as spheromaks) of CT-to- 
background magnetic field (~ 0.1). Simulations with higher injection velocity were 
made possible by the shock-handling capability of our three-dimensional (3D) ideal 
magnetohydrodynamic (MHD) code [31]. The lower CT-to-background magnetic field 
ratios, compared to past work, was enabled by the higher resolution of our code which 
allowed the boundary layer between the CT and background plasma to be properly 
resolved. 

The paper is organized as follows. In Sec. O we describe the the problem setup 
including initialization of the CT and background slab plasma and the numerical model. 
We present the simulation results in Sec. [3l and our conclusions and implications for 
future CT fueling experiments are given in Sec. HI 
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2. Problem setup and numerical model 

A low j3 CT with spherical radius = 1, centered initially at Xb = 0, i/b = 
and Zb = Zbfi = —12, is injected along the z axis into a lower density background 
plasma with injection velocity finj (see Figure [1]). The basic model assumptions and 
numerical treatments are briefly summarized here; they are essentially the same as 
those in Li et al.[32] where more details are given. This code uses high-order Godunov- 
type flnite-volume numerical methods. These methods conservatively update the zone- 
averaged fluid and magnetic fleld quantities based on estimated advective fluxes of mass, 
momentum, energy, and magnetic fleld at the zone interface [32]. The divergence- free 
condition of the magnetic fleld is ensured by a constrained transport (CT) scheme 
[33]. All simulations were performed on the parallel Linux clusters at Los Alamos 
National Laboratory. The nonlinear system of time- dependent ideal MHD equations in 
3D Cartesian coordinates {x, y, z) is given here: 

§^ + v-(p^l =0, (1) 

^^^""^ ' V- (pvv + {p+^)I-Bb] = , (2) 



dt 
dE _ 
+ ^ 



E + P+ — ] v-B{v-B) 



0, (3) 



— -Vx(t;x5) =0, (4) 

in which p, p, -u, B and E are the density, (gas) pressure, flow velocity, magnetic 
field, and total energy, respectively. I is the unit diagonal tensor. The total energy 
is = p/(7 — 1) + pf^/2 + -B^/2, where 7 is the ratio of the specific heats. Note that 
a factor of -y/ivr has been absorbed into the scaling for both the magnetic field B and 
current density j. It should be noted that the details of effects such as reconnection and 
heat evolution could not be addressed accurately due to the ideal MHD model and the 
use of a simplified energy equation. 

It is well established empirically in coaxial gun spheromak experiments that, under 
proper conditions, a spheromak "magnetic bubble" (a low (3 CT) will be formed by 
the gun discharge [M]. In our simulations, we do not model the CT formation process 
and instead start with a pre-formed CT moving toward the background slab plasma at 
speed finj. The stationary background plasma is composed of a slab plasma confined by 
a uniform background magnetic field Bp in the x direction. For simplicity, we assume 
the background plasma has uniform initial number density pp = 0.1 and uniform initial 
temperature Tp = 0.1. 

In the simulations reported here, the CT structure is similar to the one given in 
Liu et al. [35]. The density profile of the CT plasma with radius rf, = 1 is given by 

pb oc r^exp[-r^ - {z^ - Zbf], 



up to a normalization coefficient and a uniform temperature Tb, where Tc = \/ x'^ + y'^ 
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and Zc = z (see Figured]). The density profile used here has its peak shifted from the 
center of the CT, approximating a spheromak. 

The CT magnetic field is determined by three key quantities: the length scale of 
the bubble magnetic field = 2, the amount of poloidal flux and the index a, 
which is the ratio of the CT toroidal to poloidal magnetic fields. For simplicity, the CT 
magnetic field Bf, is also assumed to be axisymmetric. The poloidal flux function "ifp is 
specified as 

cx exp[-rl - {z, - z^f] . (5) 

The poloidal fields, up to a normalization coefficient, are 

1 d'^ 1 9\1> 

Bb,r^ = 56,zc = (6) 

Tc OZc Tc OTc 

while the toroidal magnetic field is 

Bb,^, = = arcexp[-r^ - {z^ - Zbf] . (7) 
rc 

The azimuthal component of the CT Lorentz force is zero, but the total azimuthal 
Lorentz force due to the combined fields and currents of the CT may be non-zero. 

The CT also has uniform injection velocity vi^j and uniform rotation angular speed 
uj. In this paper the ratio of the CT's toroidal to poloidal magnetic fields a, the 
rotation speed of the CT ui, and the specific heat 7 are taken to be a/IO, and 5/3, 
respectively. Physical quantities are normalized by the characteristic system length 
scale -Ro = 10 cm, density po = 7.77 x 10~^ g/cm~^ (corresponding to plasma number 
density, which is the sum of electron and ion number density, of 18.6 x 10^^ cm~^), and 
velocity Vq = 1.7 x 10^ cm s~^. Other quantities are normalized as: time t = 1 gives 
Ro/Vo = 5.9 X 10"^ s, magnetic field B = 1 gives {ATrpoV^Y^'^ = 5.3 x 10^ G, and energy 
E = 1 gives poV^Rl = 2.24 x 10" ergs. 

The boundary conditions are all perfectly conducting in the y and z directions 
except at the port where the CT is injected, while in the x direction outflow boundary 
conditions are employed in order to mimic the toroidal geometry of a tokamak. Since 
the toroidal dimension is much larger than the poloidal dimension in a real tokamak 
and we focus on the early stage of the CT evolution before one toroidal propagation 
time of the Alfven wave induced by the CT propagation, the out-flowing boundary 
condition is more appropriate than the periodic boundary condition used in Suzuki 
et al. [7]. In order to minimize the influence of the entrance port, the port will be 
switched on when the top of the CT reaches the bottom boundary at t = and 
it will be switched off at t = 5/^'inj after the CT has fully entered the computation 
domain {t ~ A/vi^j)- Suzuki et al. have pointed out that the boundary condition in the 
background magnetic field (x) direction is important, i.e., magnetic reconnection has 
more influence on CT deceleration with perfectly conducting boundary conditions than 
with stress-free boundary conditions such as outflow and periodic boundary conditions 
J7\. The total computational domain is |a;| < 9, \y\ < 9, and 1^1 < 9, corresponding 
to a (180 cm)^ box in actual length units (assuming the physical dimension of the 
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injected CT radius = 10 cm). The numerical resolution used here is 400 x 400 x 800, 
where the grid points are assigned uniformly in the x, y, and z directions. A cell 5x 
(= Sy = 2Sz = 0.045) corresponds to 0.45 cm. Since the plasma skin depth and ion 
gyroradius based on the parameters of ITER are no more than 6x, the simulations based 
on an MHD model are appropriate. 

3. Results 

In this section we present ideal MHD simulation results on the injection of a low /3 
CT into a hot strongly magnetized plasma. We organize our results into three topics: 
(1) parameter regimes, in terms of finj and ratio of CT-to-background magnetic field, 
of CT evolution including the identification of a promising regime for ITER-relevant 
precise core fueling, (2) detailed description of the CT evolution for the ITER-relevant 
regime, and (3) dependence of the CT penetration depth on fjnj and the initial CT 
density. 

3.1. Parameter regimes of CT injection 

Based on our simulation results, we find that the evolution of the injected CT 
depends predominantly on the initial injection speed Winj and the initial ratio of CT-to- 
background magnetic field. As shown in Figure Ej there are several qualitatively distinct 
regimes of CT injection in terms of the above two parameters. 

First, below a threshold injection speed Vl, the CT is unable to penetrate the 
background plasma at all (see left panel of Figure[3]). The conducting sphere (CS) model 
[1] requires that the initial CT kinematic energy exceeds the background magnetic field 
energy excluded by the CT volume. 



Our results show that Vl, which is inferred (~ 0.5) from extrapolation of the data 
shown in Fig. [9](a), is less than Vac- This is because of the compression of the CT 
during the penetration, which is ignored by the CS model (see discussion of §3.21) . i.e., 
the CT actually excludes less volume than its initial volume, therefore leading to a 
smaller initial injection speed threshold. And since this compression is related to the 
background plasma magnetic field (see discussion in §3.21) . this injection speed threshold 
is thus dependent on the background plasma magnetic field. However, Vac derived from 
the CS model still gives a rough estimate of the lower limit of the injection speed needed 
for penetration. 

Second, above a different threshold injection speed Vap, a strong shock and 
wavefront are observed to develop ahead of the CT, dominating the system evolution 
(see middle panel of Figure [3]). This threshold is determined by the condition when vi^j 
exceeds the Alfven speed of the background plasma, i.e., vi^j > Vap = B^j ^Jp^. This 
shock/ wavefront-dominated regime is probably not favorable for CT fueling because it 




f inj > Vac 



(8) 
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may strongly perturb the tokamak equilibrium. Also, in this case the CT deposition is 
highly nonlocal {middle panel of Figure [3]). For comparison, the result (finj = 1.1) 
without shock/ wavefront development is also presented {right panel of Figure [3]) 
(Vap = 3.16 in this case). 

Third, in the regime Vi, < finj ^ Vap, CT evolution is further determined by the 
ratio Sj. = B^/Bp of the CT field Bf, to the background plasma magnetic field Bp. If 
Sj. ^ 10, a large non-zero initial Lorentz force results in strong CT expansion that 
also leads to the development of a shock and a wavefront that dominate the system 
evolution [35]. As stated above, this does not favor controlled plasma fueling. In the 
regime 1 < S'^ < 10, which is the regime Suzuki et al. have discussed extensively [HlinilTj, 
the CT is decelerated by both the magnetic pressure and magnetic tension forces. In 
this regime, the CT tilts while reconnection occurs between the CT and background 
plasma magnetic fields. Suzuki et al. proposed the Non-slipping Conducting Sphere 
(NS) model [5],[6l[7], which matches their simulation results of vertical injection without 
a magnetic field gradient pretty well [7j. Our simulation results in this regime verify 
their conclusions. Because it is difficult for the CT field Bj, generated by a co-axial 
gun to be larger than the tokamak field Bp of several Tesla or more, it is important to 
explore the regime Sr < 1, which is the primary focus of the remainder of this paper. We 
have identified this regime ( right panel of Figure [3]) as a promising one for CT fueling 
of ITER-relevant plasmas due to the precise spatial deposition of the CT and the deep 
penetration that can be achieved for core fueling applications. In this paper, we have 
focused on background plasmas with (3p = 0.02 like in ITER, while simulation results 
with f3p = 0.2 like in the National Spherical Torus Experiment (NSTX) will be studied 
in follow-on work. 

3.2. CT Evolution 

In this sub-section, we describe in detail the CT evolution for the ITER-relevant regime 
of Vl < finj ^ Vap and Sr < 1. In this regime, the CT evolution can be divided 
into three stages: (1) initial CT penetration, (2) CT compression in the propagation 
direction {z) and reconnection, and (3) coming to rest and spreading in the x (toroidal) 
direction. Magnetic reconnection, starting late in the second stage in our simulations, 
arises due to numerical diffusion and leads to mixing between the CT and background 
plasmas. 

In order to mimic CT injection with ITER-relevance, we adopt physical quantities 
as given in Table[I]and choose the injection speed finj such that Vac = 1-0 < fmj = 1-1 < 
Vap = 3.16 (this falls into the "ITER-relevant" case of Figure [2]). Figure H] displays the 
time evolution of the plasma density (color contours in common logarithmic scale) in 
the x-z plane with y = 0. The white solid contour lines indicate the magnetic pressure 
Pb = B'^/2. Fig. [5] displays the axial profile of x and y integrated density f^f pdxdy 
at different times corresponding to Fig. HI The entire evolution can be described by the 
three stages mentioned above. 
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Figure E] displays the magnetic field B^^, (arrows) and current density jy (color 
contours) in the x-z plane at y = 0. During the first stage (initial penetration), the CT 
experiences a very strong magnetic obstacle, and therefore the CT plasma is highly 
compressed and the plasma density increases at the interface between the CT and 
background plasmas. A large plasma current also appears at the interface due to the 
compression of the background magnetic field seen in Figure [6]^ a). This current sheet 
is bent and broken into two parts (Figure [6]^b)) due to the magnetic field configuration 
of the CT field. Some reconnection takes place at the left part of the current sheet as 
shown in Figure [71 The CT is successively decelerated by the magnetic tension force of 
the background magnetic field. 

After the CT has fully entered the background plasma region (after t ~ 4), the 
background field lines, some of which have reconnected with the CT magnetic field at the 
CT's leading edge, are reconnected again with the CT magnetic field at the CT's trailing 
edge, as seen in Figure [6](c). As pointed out by Suzuki et al.[5], via this process the CT 
is separated from the background fields lines, through which the magnetic tension force 
decelerating the CT is relaxed. At t = 4.375, a magnetic configuration schematically 
shown in Figure[7]is formed. This interpretation is supported by Figure[Hl which displays 
the axial distributions of density and B^. From Figure [H we can see that B^ changes 
sign from positive to negative and then to positive again. The transition from positive 
to negative happens at larger z on the left hand side while at smaller z on the right 
hand side, which is the case in FigmeUlright) (see two dashed lines in Figure [7](h^/it)). 
From Figure [TJ primary reconnection sites are at the upper left and lower right sections 
of the CT. The reconnection process allows the high-density CT plasma to escape from 
the CT and eventually flow outward along the background magnetic field horizontally 
(Figure [7]). The CT plasma starts to contract in z and expand in x. The reconnection 
is asymmetric about the CT axis, and this asymmetry results in CT plasma outflow 
in a direction that is not completely in the x direction, but rather obliquely (panel (d) 
of Figure [6]). Eventually it will become more parallel to x {right panel of Figure [3]). 
Magnetic flux is being destroyed as well. During this process the MHD wave drag from 
Alfven waves induced by horizontal plasma outflow lU H] might further slow down the 
CT plasma. Note that Suzuki et al.[7j showed that the CT penetration depth, based 
on a model with magnetic tension force as the main deceleration mechanism, matches 
simulation results very well, implying that MHD wave drag forces may not be important 
in CT deceleration. The initially injected magnetic and perpendicular kinematic energies 
are converted into parallel kinematic energy. Contrary to [Ij and O [6l [7] , CT tilting, 
the time scale of which is proportional to ^y BbBp/ pb, is not observed in our case due to 
the fast injection and short CT transit time as required for the ITER-relevant regime 
(Table [ID . 

After the high-density CT plasma has been depleted during the compression stage 
(after t ~ 9), the CT field diminishes and only the slightly perturbed background field 
survives, almost orienting to the direction of the initial background magnetic field (see 
Figure [Hand Figure [6]^d)). Almost all of the initially injected magnetic and kinematic 
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energy now survive in the form of parallel kinematic energy. The CT comes to rest in 
the injection direction {z). A narrow elongated structure along x results, as seen in 
the lower middle and lower right panels of Figure HI This line-shaped structure with a 
spread of only Az ~ 0.2 (much smaller than the size of the background plasma, see also 
Fig. [5]) implies a very precise fuel deposition capability. 

3.3. CT penetration depth 

In this sub-section, we establish the dependence of CT penetration depth, an important 
parameter for CT fueling of tokamaks, on the experimentally controllable parameters 
of CT injection speed finj and initial density pb, which collectively determine the initial 
CT energy. The penetration depth S is defined as the axial {z) distance between the 
final mean position of the injected CT and the injection location on the boundary. 
We find that S is highly dependent on the initial injection speed finj and CT density 
Pb- Figure [9](a) displays the relationship between the S and finj, showing that S is 
proportional to finj if Vl < finj < Vap. Figure [9]^b) displays the relationship between the 
CT density and penetration depth, which shows that the penetration depth S increases 
with the CT density. The CT spread in the z direction is around ~ 0.2 in all cases (see 
also Fig. [5]), which is very small compared to the background plasma size ~ 20. These 
two empirical relationships are important in the sense that they provide clues for how 
to choose the injection speed and CT density to get precise fuel deposition, therefore 
controlling the core plasma profile in a large tokamak such as ITER. It is very hard to 
manage that by other methods such as pellet injection |26j . 

For the case shown in Figure Mjo), the parameters given in Table [H are adopted 
(winj = 2.0) except CT density pb and CT field Bb. We choose them such that the 
plasma parameter f3 = 0.2 is kept constant for all the simulation results shown in 
Figure [9](b). We keep (3 constant since experimentally this is more reasonable than 
changing CT density or CT field strength independently. Thus, increasing CT density 
means increasing both the initial kinematic energy Ek = J l/2pbv'^^^dV , where dV is 
the infinitesimal CT volume, and initial magnetic energy E^ = J l/2B^dV. However, 
because of the limiting criterion given in Eq. [HI the initial CT kinematic energy Ek 
should be much larger than the initial CT magnetic energy Em since the background 
field is much larger than the CT field in all simulation results presented in Figure Mjo) 
{Bb/Bp = 0.1). Also the residual initial total force j x B — Vp is proportional to 
the CT density ~ 0{pb). Therefore, increasing the CT density pb would elevate the 
initial residual force, which would result in shock/wavefronts. Thus, with too large a 
density ratio {pb/pp ^ 15, which is inferred from simulation results), the fueling is not 
as localized as the case with smaller CT density, while generating similar phenomena as 
the shock/ wavefront dominated cases. 

Results of Figure [9] remind us of the importance of initially injected energy upon 
the penetration depth. As discussed in §3.21 how long the compression stage lasts is 
determined by the sum of the initially injected kinematic and magnetic energy. The 
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larger the sum is, the longer it will take for the reconnection to dissipate/convert the 
energy. If we assume that the deceleration mechanism of the CT is mainly due to the 
magnetic tension force, which is independent of the initially injected energy [7], then 
the larger initially injected energy would lead to deeper fueling given the same injection 
speed finj (since plasma parameter -C 1, the internal energy of the CT plasma can be 
ignored). Certainly the larger injection speed would lead to deeper penetration. 

Figure [10] shows the time evolution of the net toroidal magnetic flux tpt = J BydS 
(only positive By is selected) with t>inj = 1.1 and ph/ Pp = 10. Before t = 6.25, the net CT 
toroidal magnetic flux increases due to the initial penetration and compression of the 
magnetic field lines at the interface of the CT and background plasma. After t = 11.25, 
the line structure reaches the horizontal boundaries and some toroidal magnetic flux 
flows out. Between t = 7.5 and t = 11.25, the only source of the destruction of the net 
toroidal magnetic flux is magnetic reconnection due to the numerical diffusion. If we 
fit the decay as iptif) / iptit = 0) = exp(— t/rres), the "resistive" dissipation time due to 
numerical diffusion is Tres ~ 10.9. Therefore the global resistive decay due to numerical 
diffusion is not important on the time scales of CT transit time [t ~ 5) given injection 
speed Winj = 1.1. In the simulation, mixing arises from numerical diffusion induced 
reconnection and occurs on the same time scale as the CT transit time. However, the 
mixing should be less of a factor in reality due to the much smaller diffusion in a high 
temperature tokamak. 

4. Conclusions & discussion 

In this paper we presented nonlinear ideal MHD simulation results of a compact toroid 
injected into a hot strongly magnetized plasma. The simulations are intended to provide 
insights into CT fueling of an ITER-class tokamak. As a first step, we have investigated 
the problem of a high density low (3 CT injected into a slab background plasma with 
uniform magnetic field and density. We intend to investigate the injection of high j3 
CT's as well as unmagnetized dense plasma jets, and also incorporate more realistic 
background profiles in follow-on work. 

Our main findings are as follows. A regime is identified in terms of CT injection 
speed and CT-to-background magnetic field ratio that appears promising for precise core 
fueling. Shock-dominated regimes, which are probably unfavorable for tokamak fueling, 
are also identified. The CT penetration depth is proportional to the CT injection speed 
and density. For the regime identified as favorable for precise core fueling, the entire CT 
evolution can be divided into three stages: (1) initial penetration, (2) compression in the 
direction of propagation, and reconnection with the background plasma, and (3) coming 
to rest and spreading in the direction perpendicular to injection. Tilting of the CT is 
not observed due to the fast transit time of the CT across the background plasma. 
Reconnection occurring at the upper left and lower right portions of the CT fragments 
the CT and leads to CT plasma outfiow horizontally along the background magnetic 
field lines, forming a line-shaped structure almost parallel to the background magnetic 
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field. We have also studied how the penetration depth is determined and have found an 
empirical relationship between the penetration depth versus the initial injection speed 
and CT plasma density. The penetration depth is proportional to the injection speed 
and CT density. When the CT comes to rest, it is highly compressed in the direction 
of injection, implying that very precise deposition is possible. 

In a device with high bootstrap current fraction, optimized density and pressure 
profiles must be maintained. A core fueling system is all that a burning plasma may 
have to alter core plasma conditions and for burn control. Initial density peaking via 
core fueling provides more flexibility to reach ignition. If a CT injection system is within 
the ITER-relevant, no-shock/wavefront regime as required by finj <^ Vap, Pb/Pp ^ 15 
and Bh/Bp <^ 1, then precise core deposition with minimal background equilibrium 
perturbation may be possible. 

Because Vac ~ 1700 km s~^ for ITER (using parameters in Table [1]), it will be a 
challenge for present CT injectors with f ~ 400 km s~^ to achieve deep penetration. A 
potential solution is to increase pb/pp to at least 100, which would reduce Vac to a more 
reasonable value (~ 170 km s~^). Due to computational limitations, a simulation of such 
a case (very dense CT) is beyond the scope of this study. However, to provide insight 
into the penetration characteristics of a very high density injected plasma, we report here 
some preliminary results on the injection of a dense {pb/Pp = 538 or ~ 10^^ cm~^), 
cool {Tb/Tp = 3.33 X lO"'' or = 2.5 ev), and very weakly magnetized plasma [37] , 
which would be an interesting alternative to CT injection. Figure [TT] suggests that this 
type of injected plasma is capable of reaching the core of the tokamak, albeit with worse 
localization as discussed in Sec. 13. 3[ This will be reported on in more detail in follow-on 
work. 
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Figure 1. Schematic of the simulation geometry showing the coordinate system. In 
the text, the direction along the z— axis is defined as the axial direction. 
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Figure 2. Qualitative behavior for CT injection in terms of CT injection speed and 
magnetic field strength. 
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Figure 3. (color) Density in the x-z plane at y ~ 0. Left panel: at i = 20 with 
injection speed Wi^j = 0.3, which shows no penetration. Middle panel: at t = 2.5 
with injection speed Winj — fO.O, which shows non-localized deposition. Right panel: 
at t = 20 with injection speed Winj = f.f, which shows highly localized deposition. 
Vav = 3.f6. 



Table 1. Normalized physical quantities. Note that Vac = Bp/ ^Jpl = 1.0. Injection 
speed Winj is between and Vap . Plasma is assumed to be composed of half deuterium 
and half tritium. The plasma density is the sum of electron and ion densities. 
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numerical 


physical 


numerical physical 


Magnetic Field 


5, = 0.1 


0.53 T 


Bp = 1.0 5.3 T 


Density 


Pfe = 1.0 


1.86 X 10^^ cm-3 


Pp = 0.1 1.86 X 10^^ cm-3 


Temperature 


Tb = 0.001 


75 ev 


Tp = 0.1 7.5 kev 


plasma (3 = 2pT / < > 


&b 


= 0.2 


Pp = 0.02 


Alfven Speed Va = Bj ^ 


^AB = 0.1 


170 km s"^ 


Vap = 3.16 5.4 x 10^ km s"^ 


Sound Speed K = 


Kb = 1.7 X 10-^ 


2.89 km s^i 


Kp = 0.41 7.0 X 10^ km 
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Figure 4. (color) Density (logarithmic scale) in the x-z plane at y = as a function 
of time. The white lines are contours of magnetic pressure pb = B'^ /2. For t = 2.5, 
PB e [0.0, 1.6843]; for t = 3.75, pb e [0.0, 1.7612]; for t ^ 7.5, ps G [0.0, 1.4405]; for 
t = 8.75, pB e [0.0, 1.1872]; for t = 12.5, pb e [0.0, 1.0219]; for t = 15, pb G [0.0, 1.008]. 
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Figure 5. (color) x and y integrated density pdxdy versus z at different times 
corresponding to Fig. 01 Black dash: t = 2.5, black dash dot: t = 3.75; black solid: 
t = 7.5: red dash: t = 8.75; red dash dot: t = 12.5, red solid: t = 15. 
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Figure 6. (color) Magnetic field and current density jy in the x-z plane as a function 
of time at ?/ = 0. The color contours show jy while arrows show and -Bz, which 
are normalized to their maximum values, (a) t = 2.5; (b) t = 4.375; (c) t = 6.25; (d) 
t = 8.75. Please note that the scales for the abscissa and ordinate are not identical. 
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Figure 7. (color) Diagram of the magnetic configm-ation in the contraction stage. 
Red color indicates the high-density CT plasma. Yellow lines indicate magnetic field 
lines, while blue arrows indicate CT plasma flow patterns. See Figure [8] for dotted 
black lines. 
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Figure 8. Axial profiles of several quantities at t = 4.375 with {x,y) — (—1.125,0) 
{left) and {x,y) = (0.855,0) (right) . The density p and magnetic field strength 
in the x-direction are shown for evaluating the magnetic configuration shown 
in Figure fright). These two axial profiles correspond to the two dash lines in 
Figure [7^ri(7/it). 
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Figure 9. (a) Penetration depth S versus injection speed Vi^^ with pb/ Pp = 10 and 
(b) Penetration depth S versus CT density ph with Ui„j = 2.0. Both panels assume 
Bb/Bp = 0.1. The penetration depth is defined as the distance between the final mean 
position of the injected plasma and the bottom boundary. The error bar indicates the 
final size of the injected CT in the ^-direction. In panel (a), the dash line indicates 
Vl ^ 0.5 and the dash dot line indicates Vap = 3.16. 
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Figure 10. Decay of the net toroidal magnetic flux ipt = J BydS, where only 
positive By is selected with Vinj = 1-1 and Pb/Pp = 10. The dashed Hne fits the 
data between t = 7.5 and t = 11.25 with the formula ^t{t)/'^t{'t = 0) = exp {—t/Tj-es), 
where Tres = 10.9 is the resistive dissipation time due to numerical diffusion, and it is 
much larger than the CT transit time 5 given injection speed Vinj = 1.1. 
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Figure 11. (color) Density in the x-z plane of y = at t = 150 with injection speed 
"^inj — 0.12, pb/pp = 538, and Bt/Bp = 0.001, which shows worse localized deposition. 
Please note that the CT is still moving with almost the initial injection speed. 



